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I.  INTRODUCTION 

Visual  Flight  Simulators 

Visual  Flight  Simulators  are  used  to  enhance  the  training  of  pilots  by 
adding  visual  cues  to  the  simulated  mechanical  and  physical  aspects  of  a  par¬ 
ticular  aircraft.  The  ultimate  visual  goal  is  to  provide  a  realistic  view  of 
an  environment  about  the  simulated  aircraft  to  Increase  training  effective¬ 
ness.1  Imagery  is  often  provided  by  Computer  Image  Generation  (CIG)  and  dis¬ 
played  via  video  projection  systems.  Typically  the  CIG  system  takes  account  of 
viewpoint  and  heading  direction  within  a  mathematically  modeled  landscape 
(database)  to  generate  a  view  of  this  database  during  a  simulated  flight.2 
The  imagery  generated  is  in  a  video  format,  usually  a  1023  line  raster,  and 
can  be  projected  by  light  valves  on  the  screen  surface.3  Image  content  depends 
on  the  database  model  and  may  Include  a  horizon  line.  A  trainee's  view  of  this 
imagery,  in  conjunction  with  the  aircraft  dynamics  and  structure,  can  induce 
many  physical /psychological  effects  of  actual  flight  depending  on  the  attention 
paid  to  modeling  the  actual  aircraft  and  database.4 

Wide-Angle  Visual  Displays 

The  goal  of  realism  in  visual  displays  has  led  to  the  use  of  very  wide 
angle  displays  filling  a  horizontal  field  of  view  (FOV)  of  90°  or  more  at  the 
pilot's  viewpoint.5  This  has  led  to  the  use  of  spherical  screens  (domes)  with 
a  number  of  projectors  filling  parts  of  the  FOV  to  form  a  wide-angle  scene.6 
Ideally,  the  exit  pupils  of  the  projectors,  along  with  the  viewpoint,  should  be 
located  at  the  center  of  the  spherical  screen  to  eliminate  distortions.  How¬ 
ever,  due  to  the  physical  structure  of  the  simulated  aircraft  and  the  size  of 
the  projection  systems,  the  exit  pupils  cannot  all  be  placed  at  the  screen 
center.  Also,  the  viewpoint  and  projector  exit  pupils  cannot  coincide.  Neces¬ 
sarily,  this  forces  the  oblique  projection  of  the  Imagery  onto  the  dome  result¬ 
ing  in  distortion  of  the  Imagery. 
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Distortion 


Distortion  refers  to  the  geometry  of  ar,  image  as  compared  to  the  geometry 
of  the  objects  involved.  Perhaps  the  best  way  tc  describe  the  concept  of  dis¬ 
tortion  as  related  to  optical  systems  is  to  consider  first  the  concept  of  a 
distortionless  lens  mapping.  Figure  1  shows  the  mapping  process  for  a  well 
corrected  F  Tan  e  lens  as  a  rectangular  object  is  mapped  through  the  lens. 

The  resulting  image  formed  on  the  image  plane  is  again  a  rectangular  figure. 

The  radial  distance,  R,  to  an  image  point  cn  this  image  plane  is  determined  by 
the  tangent  of  the  angle  that  an  object  point  subtends  from  the  optical  axis 
of  the  lens  system.  Hence,  the  lens  mapping  equation,  R  =  F  Tan  c,  where  F  is 
the  focal  length  of  the  lens.  If  the  direction  cf  the  mapping  process  is  re¬ 
versed,  we  now  have  the  case  cf  the  lens  being  used  as  a  projection  lens  in¬ 
stead  of  use  as  a  taking  lens.  For  this  case,  any  imagery  placed  on  the  image 
plane  will  be  transferred  to  the  object  plane  without  distortion.  This  is  the 
concept  of  a  distortionless  lens  mapping. 

Next,  we  should  consider  the  concept  of  distortionless  viewing  of  a  pro¬ 
jected  image.  If  it  were  physically  possible  to  place  the  eye  of  a  viewer  at 
the  exit  pupil  of  a  F  Tan  c  lens  as  in  Figure  2,  then  the  viewer  would  perceive 
r.c  distortion  of  the  projected  imagery.  However,  if  the  viewer's  e>e  is  re¬ 
moved  from  the  exit  pupil  location,  the  shape  of  a  rectangular  figure  would  nc 
longer  appear  rectangular.  This  removal  of  the  eyepcint  from  the  lens  e>it 
pupil  results  in  a  form  of  distortion  known  as  Perspective  Distortion,  which 
usually  results  ir.  rectangular  objects  exhibiting  e  keystone  shape.  The  key 
zo  reouc.ng  Perspective  Distortion  is  tc  place  the  exit  pupil  of  the  projector 
as  close  as  possible  tc  the  viewer's  eyepoint. 

In  a  dome  display  system  there  exists  another  form,  cf  distortion.  This 
form  is  labeled  as  Geometric  Distortion,  which  involves  the  projection  of 
imagery  onto  curved  screen  surfaces  at  oblique  angles.  The  result  is  that  the 
projection  of  straight  lines  onto  the  screen  surface  are  viewed  as  being 
curved  lines  by  the  viewer.  The  amount  of  distortion  produced  depends  on  the 
size  end  shape  of  the  screen,  and  projector/viewpoint  positioning.  As  the 
viewpoint  and/or  the  projector  are  displaced  from  each  other,  or  removed  frer 
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Figure  2.  Distortionless  viewing 
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the  center  of  the  screen,  the  apparent  size  and  shape,  as  well  as  the  angular 
subtense  and  position  of  objects  will  vary  accordingly.  It  should  be  noted 
that  if  the  viewpoint  and  lens  exit  pupil  are  coincident,  then  the  shape  of 
the  screen  cannct  contribute  to  distortion7.  It  is  the  purpose  of  the  com¬ 
puter  program  to  be  described  to  consider  these  two  forms  of  distortion.  Per¬ 
spective  and  Geometric  Distortion,  in  spherical  screen  displays. 

Distortion  definitions  vary  according  to  the  type  cf  projection  system 
involved.  The  method  of  calculating  distortion  in  this  report  is  based  on  the 
Institute  cf  Radio  Engineers  (IRE)  Standards  for  Television.8  This  method  de¬ 
fines  the  Geometric  Position  Error  (GPE)  for  any  point  on  a  flat  display 
screen  as  the  magnitude  of  the  distance  from  the  point  to  its  ideal  location. 
This  implies  a  radial  distance  measurement  front  the  ideal  location  of  a  point 
to  its  actual  position.  The  percentage  of  distortion  is  then  found  by  divid¬ 
ing  the  GPE  by  the  full  field  height  of  the  image. 

The  use  of  video  projection  systems  offers  the  option  of  correcting  dis¬ 
tortions  by  altering  the  scanning  geometry  of  the  video  imaging  system  (i.e., 
raster  shaping).  In  this  way,  the  image  can  be  pre-distorted  (as  compared  to 
a  rectangular  raster)  before  projection,  according  to  specific  projector/ 
viewpoint  relationships,  in  order  to  produce  a  non-distorted  view.9  Alter¬ 
nately,  the  object  points  can  be  remapped  in  the  CIG  computation  plane  to 
provide  the  required  object  pre-distortion  before  being  placed  on  the  rectangu¬ 
lar  raster  for  projection.  Raster  shaping  and  remapping  nay  be  combined,  if 
desired,  to  reduce  the  complexity  of  the  individual  corrections.  These  correc¬ 
tions  must  take  into  account  the  type  of  projection  lens  to  be  used  and  Its 
inherent  distortion,  for  the  lens  is  an  integral  part  of  the  system. 

Lens  Mappings 

In  the  computer  program  that  will  be  described  there  are  four  types  of 
lens  mappings  considered  for  the  projection  lens.  They  are: 

1.  F  Tan  ©  (distortionless  lens), 

2.  F  T^n  0  with  primary  distortion. 


5 


NAVTRAEQU I PCEN  IH-332 


3 .  F  £, 

4.  F  Sin  o. 

All  these  lens  mappings  imply  a  radially  symmetric  mapping,  with  the 
center  of  the  mapping  plane  on  the  optical  axis  of  the  lers  system.  The 
F  Tan  0  lens  provides  a  one-to-one  mapping  fron  object  space  tc  image  space 
with  a  magnification  constant  included.  The  F  Tar.  0  lens  places  an  irage 
point  cn  the  image  plane  according  tc  the  tangent  of  the  angle  between  the 
optical  axis  and  the  object  point.  The  F  Tan  c  with  Primary  Distortion  is 
defined  as  a  departure  from  F  Tan  t  due  to  an  approximation  of  Tan  t  by  only 
two  terms  of  a  power  series,  the  resulting  equation  is: 

Rp  =  Rt  (1  +  ( DFACTGR  *  Rt2)), 

where  R  is  the  radial  aistance  for  the  primary  distortion  mapping,  R.  is  the 
P  ^ 

radial  distance  for  Tan  0  mapping  and  CFACTOR  is  the  primary  distortion  fac¬ 
tor.  An  F  0  lens  maps  object  space  tc  image  space  according  to  the  angle  ir. 
radians  to  an  object  point  from  the  optical  axis,  while  the  F  Sin  0  implys  a 
mapping  according  to  the  sine  of  the  angle  from  the  optical  axis.  Appendix  B 
contains  the  derivation  for  the  F  Tan  &  mappings  and  the  conversion  to  the 
three  other  mappings  mentioned  above. 

For  the  points  on  these  mappings,  the  angle  from  the  optical  axis  of  the 
projection  lens  can  be  found  by  calculating  the  radial  distance  to  the  point 
in  question  and  then  converting  the  radial  distance  tc  an  angle  depending  on 
the  type  of  mapping. 

Figure  3  shows  the  mapping  process  for  a  F  Tan  g  lens  frcm  the  View  Window 
on  the  dome,  through  the  View  Plane,  and  onto  the  Target  Plane  of  the  ler.s. 

The  figure  also  depicts  the  viewpoint  mapping  of  the  View  Window  onto  the  Com¬ 
parison  Plane.  This  mapping  produces  a  two-dimensional  representation  of  the 
View  Window  as  seen  from  the  viewpoint. 

The  projection  mappings  produced  by  the  program  are  mappings  at  the  View 
Plane,  which  provide  a  perspective  view  of  the  object  points  on  the  dome  from 
the  projection  point.  A  perspective  view  cf  the  Target  Plane  frcm  the 
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Figure  3.  Projection/Viewing  Diagram 
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projection  point  (a  view  of  the  required  raster  pattern)  is  obtained  by  a 
rotation  of  the  mapping  about  its  center  by  180°  and  observing  the  pattern 
through  the  reverse  side  of  the  mapping.  In  all  mappings,  the  value  of  F  is 
considered  to  be  unity,  allowing  mappings  for  a  lens  of  focal  length  f1  to  be 
represented  by  multiplying  the  mapping  point  coordinates  by  f1. 

The  Comparison  Plane  is  a  perspective  view  of  the  View  Window  from  the 
viewpoint  and  contains  the  image  of  the  non-distorted  reference  raster  which 
has  been  placed  on  the  View  Window.  Use  of  the  Comparison  Plane  coordinates 
requires  the  scaling  of  the  coordinates  to  reduce  the  Comparison  Plane's 
height  to  the  height  of  the  View  Plane.  The  scaling  factor  is  defined  as  the 
ratio  of  the  View  Plane  height  to  the  Comparison  Plane  height.  For  purposes 
of  distortion  calculations  in  the  program  to  follow,  the  GPE  is  found  by  ef¬ 
fectively  overlaying  the  View  Flane  on  the  Comparison  Plane  via  coordinate 
tatles  for  the  planes.  The  distance  from  a  point  on  the  Comparison  Plane  (XC, 
YC)  to  a  point  on  the  View  Plane  (X,  Y)  is  found  via  the  equation: 

D  =  [(X  -  XC)2  +  (Y  -  YC)2]l/2. 
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II.  THE  PROGRAM  -  MAPTAG 

Introduction 


MAPTAG  is  an  acronym  for  Mapping  Table  and  Graph,  a  program  written  in 
FGRTRAN  for  distortion  analysis.  The  facilities  for  operation  of  the  program 
are  part  of  the  NAVTRAEQU I PCEN  Computer  Systems  Laboratory.  The  computer 
system  is  a  VAX-11/780  with  graphics  provided  on  a  Tektronix  graphic  terminal. 
Model  4014-11.  The  program  is  designed  to  find  the  required  distorted  raster 
shape  for  projection  onto  a  dome  from  a  particular  projection  point.  This 
projected  raster  is  to  be  viewed  at  a  viewpoint  as  a  non-distorted  raster. 

The  location  of  the  viewpoint  and  projector  can  be  located  inside  or  outside 
the  dome  of  radius  R.  The  raster  to  be  viewed  is  created  as  an  Input  Plane 
with  variable  height  and  width,  as  required,  to  fill  a  desired  vertical  and 
horizontal  angular  field  relative  to  the  viewpoint.  This  Input  Plane  is 
subsequently  mapped  on  the  surface  of  the  dome  to  provide  the  View  Window. 

In  order  to  describe  the  projection/viewing  system,  the  location  cf  the 
prcjectcr,  viewpoint,  and  image  points  are  referenced  to  a  3-D  coordinate 
system  located  at  the  center  of  the  dome  (Figure  4).  In  this  system,  the 
Z-axis  is  positive  upwards,  the  X-axis  is  positive  forward,  and  the  Y-axis 
is  positive  to  the  left.  Additionally,  the  angles  for  projection  end  viewing 
are  spherical  angles  referenced  to  the  positive  X-axis.  Vertical  angles  ere 
positive  above  the  dome  horizon  and  negative  below  the  horizon  with  a  maximum 
magnitude  of  90°.  Horizontal  angles  are  positive  for  a  counter-clockwise 
rotation  (positive  X-axis  into  the  positive  Y-axis)  when  viewed  from  a  point 
on  the  positive  Z-axis.  The  horizontal  angles  have  a  maximum  magnitude  of 
180°.  There  are  three  other  coordinate  systems  involved  in  the  program. 

These  are  two  3-D  coordinate  systems,  located  at  the  viewpoint  and  projection 
point,  and  a  2-D  coordinate  system  used  to  describe  the  raster  plane  defined 
as  the  Input  Plane.  This  plane  is  used  to  provide  a  non-distcrted  raster 
pattern  for  comparison  purposes. 


'  / 
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View  Window 


The  View  Window  is  defined  as  that  part  of  the  available  FOV  that  is 
being  filled  by  one  projector.  Ideally,  the  viewer  will  see  objects  project¬ 
ed  onto  this  window  as  non-distorted.  If  a  crosshatch  pattern  is  the  project¬ 
ed  object,  then  the  viewer  will  want  to  see  the  crcsshatch  pattern  with  equal 
spacing  between  the  horizontal  and  vertical  lines  on  the  view'  window.  This  is 
the  Input  Plane  and  is  constructed  by  considering  the  desired  angular  height 
and  width  as  well  as  the  position  cf  the  center  of  the  view  window  (PCTR) 
relative  to  the  viewpoint.  These  values  are  used  to  find  the  height  and  width 
of  the  crosshatch  pattern  in  dimensional  units  that  are  fixed  with  respect  to 
cimersional  units  used  for  all  3-C  locations. 


The  Input  Plane  has  a  2-D  coordinate  system  with  its  origin  set  at  the 
center  of  the  plane  (Figure  5).  In  this  coordinate  system,  the  Y-axis  is 
positive  upward  and  the  X-axis  is  positive  to  the  left.  Once  the  height  and 
width  of  the  Input  Plane  are  known,  then  the  first  point  on  the  plane  is 
found  by  dividing  the  height  and  width  in  half.  This  first  point  is  defined 
as  the  top  left  point  on  the  plane.  Subsequent  points  on  the  plane  are  found 
by  considering  the  number  of  horizontal  lines  (NHORIZ)  and  the  number  of 
vertical  lines  (NVERT)  on  the  plane.  The  width  of  the  plane  is  divided  by 
(NHORIZ-1)  to  find  the  linear  increments  along  the  lines  (HINCREMENT),  while 
the  height  is  divided  by  (NVERT-1)  to  find  the  linear  increments  between  lines 
(VINCREMENT).  To  find  the  next  point  along  a  line  on  the  plane,  the  HINCRE¬ 
MENT  is  subtracted  from  the  previous  points'  X-coordinate  with  the  Y-coordi- 
nate  remaining  the  same.  At  the  end  of  a  line,  the  VINCREMENT  is  subtracted 
from  the  previous  points'  Y-coordinate,  and  the  X-ccordinate  is  reset  to  the 
X-ccordinate  of  the  first  point.  In  this  way,  lines  are  drawn  from  left  to 
right  and  top  to  bottom  on  the  Input  Plane. 

After  a  point  on  the  Input  Plane  is  found,  the  Input  Plane  is  aligned  to 
be  normal  to  the  viewers'  Line  of  Sight  (LOS)  and  translated  to  the  desired 
View  Window  center  (Figure  6).  This  operation  is  performed  by  the  sub¬ 
routine  TROPLN  (see  Appendices).  At  this  point,  the  View  Window  is  projected 
onto  the  surface  of  the  spherical  screen  but  still  appears  to  be  a  flat  non- 
distorted  raster  from  the  viewpoint. 
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Spherical  Screen  Projection 

To  find  the  intersection  of  a  ray  projected  from  the  viewpoint  towards 
the  spherical  screen,  the  subroutine  SPTFAN  "Sphere  Point  from  Angles"  is 
used.  This  subroutine  uses  spherical  angles  of  projection  to  define  the 
direction  of  individual  rays  at  the  viewpoint.  This  distance  fron  the  view¬ 
point  to  the  screen  is  calculated  and  used  to  find  the  terminus  of  the  ray  in 
3-D  coordinates. 

The  angles  to  points  on  the  View  Window  plane  from  the  viewpoint  are 
found  by  subroutine  SAN6FPT  "Spherical  Angles  from  Points"  and  then  are  input 
to  SPTFANG.  The  View  Window  is  then  mapped  onto  the  screen  surface  and  we 
are  now  ready  to  find  a  perspective  view  of  the  View  Window  from  the  projec¬ 
tion  point. 

Perspective  View 

Another  subroutine,  PSPECTIVE,  is  used  to  find  the  perspective  view  of 
the  View  Window  from  the  projection  point.  Initially,  the  subroutine  trans¬ 
lates  the  origin  of  the  3-D  coordinates  for  the  View  Window  to  the  projection 
point.  Then,  the  View  Window  is  rotated  about  the  projector  to  align  the 
center  of  the  View  Window  (PCTR)  with  the  projector's  X-axis.  Having  accom¬ 
plished  these  operations,  the  tangents  to  points  on  the  View  Window  can  be 
found  relative  to  PCTR. 

PSPECTIVE  allows  this  tangent  mapping  to  be  altered  by  a  remapping 
through  one  of  the  three  other  lens  mappings  listed  previously.  Only  the 
tangent  mappings  will  give  a  correct  perspective  from  the  projection  point, 
with  the  three  additional  mapping  choices  distorting  the  view  according  to 
their  mapping  functions.  These  mapping  coordinates  are  sent  to  plotting 
programs  which  draw  the  View  Plane  on  a  graphic  terminal  with  hard-copy 
available. 

Reference  circles,  indicating  total  Flelds-of-View  (FOV)  of  90°  and  110°, 
are  Included  on  the  output  mappings.  These  FOV  are  generated  by  subroutine 
FOV  and  are  dependent  on  the  type  of  lens  mapping  specified.  These  reference 
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circles  can  be  used  to  determine  the  location  of  object  points  relative  tc 
the  optical  axis  of  the  lens.  Points  that  lie  on  a  circle  are  defined  as 
being  located  at  the  half-angle  of  the  total  FOV.  For  example,  if  a  point  is 
on  the  90°  FOV  circle,  then  the  angle  between  the  optical  axis  and  the  ray  to 
that  point  Is  45°.  Points  outside  the  FOV  reference  circle  are  at  greater 
angles,  while  points  Inside  the  FOV  reference  angles  are  at  an  angle  less 
than  the  FOV  half-angle.  The  FOV  reference  circles  give  an  excellent  way  to 
approximate  the  total  FOV  required  for  a  particular  projection  arrangement. 

Mapping  Tables 

The  coordinates  for  the  View  Plane  output  mapping  can  also  be  placed  in 
table  form.  Each  table  contains  specific  information  on  the  projection/ 
viewing  system  and  the  type  of  mapping  used.  The  coordinates  are  printed  in 
the  table  according  to  their  relative  position  on  the  view  plane.  Also  placed 
in  the  table  output  are  the  coordinates  of  the  Comparison  Plane,  allowing  dis¬ 
tortion  calculations  to  be  made  from  the  table.  If  table  output  is  desired, 
only  12  horizontal  points  should  be  requested  on  the  crosshatch  pattern.  This 
limit  is  due  to  the  fixed  width  of  the  printer  paper  and  the  width  of  numerals 
for  each  coordinate  value  in  the  table. 

Subroutine  TBLDIST  uses  the  coordinates  sent  to  the  mapping  tables  to 
calculate  distortion  percentages  for  nine  points  in  the  field.  These  nine 
points,  assuming  an  odd  number  of  points  or.  the  plane,  are  the  four  corner 
points,  the  central  point,  and  the  four  points  at  the  midpoint  of  each  edge 
of  the  plane.  The  percentages  found  are  included  on  the  mapping  table  output. 
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III  PROGRAM  OPERATION 

Running  the  Program 

In  order  to  operate  the  program,  variables  describing  the  projection/ 
viewing  system  must  be  entered.  Upon  instructing  the  computer  to  execute  the 
program  (RUN  MAPTAG),  the  user  is  prompted  to  enter  the  necessary  values. 
Definition  of  the  view  window  starts  with  the  location  of  the  center  of  the 
view  window  (PCTR).  This  location  is  in  3-D  coordinates  relative  to  the 
screen  center  and  must  be  on  the  screen  surface.  Next,  the  height  and  width 
in  degrees  and  the  number  of  points  across  and  down  the  window  are  entered. 

The  maximum  height  or  width  of  the  window  is  less  than  160°  due  to  the  method 
of  creating  the  input  plane.  Additionally,  the  number  of  points  should  be  odd 
for  use  in  the  distortion  calculation  subroutine. 

Describing  the  rest  of  the  system  requires  entering  the  radius  of  the 
screen  (R),  the  projector  location  EX(FX,  FY,  FZ),  and  the  viewpoint,  EN(XE, 
YE,  ZE) .  The  units  for  these  coordinates  are  arbitrary,  but  must  the  same  for 
all  3-D  coordinates.  As  an  option,  the  projector  axis  can  be  offset  from  PCTR 
during  operation  of  the  program  by  entering  degree  offsets  ether  than  zero  for 
VCFFP  and  HGFFP. 

The  rest  of  the  values  to  be  entered  are  concerned  with  the  type  of  output 
mapping  and  table.  The  type  of  lens  mapping  is  chosen  by  entering  a  value 
from  1  to  4  corresponding  to  the  four  lens  mappings  available.  If  a  Ten  c-> 
mapping  with  a  primary  distortion  is  desired,  then  the  user  is  prompted  to 
enter  a  distortion  factor  (DFACTOR).  Additionally,  for  initialization  of  the 
graphics  terminal  screen,  the  user  is  instructed  to  enter  a  graphic  scaling 
factor  (FSCALE),  screen  magnification  factcr  (SNAG),  and  offsetting  values  for 
the  origin  of  the  screen  (XEIAS,  YBIAS)  in  inches.  Finally,  the  user  can 
decide  if  table  output  is  wanted  during  the  present  run  of  the  program. 

After  the  program  has  completed  the  graphics,  the  user  is  instructed  to 
enter  "C"  to  continue  execution.  The  user  is  then  offered  the  option  of  chang 
ing  the  projection  axis  offset,  graphic  scaling,  graphic  screen  magnification, 
and  table  output  option.  If  a  change  is  desired,  the  program  again  prompts 
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the  user  to  enter  values  and  the  graphic  screen  is  cleared  in  order  to  draw  a 
new  perspective  view  according  to  the  new  values.  If  no  changes  are  desired, 
the  program  will  clear  the  graphic  screen  and  execution  is  stopped. 

Sample  Projection  Systems 

As  examples,  the  results  of  two  sample  projection/viewer  systems  are  in¬ 
cluded  in  the  Appendices.  The  first  sample  is  a  simple  projector/viewer 
arrangement  where  the  viewer  is  placed  at  the  center  of  a  20-foot  radius  dome 
and  the  projector  is  located  at  the  3-D  location  (0.0,  0.P,  12.0)  inches. 

This  results  in  a  12-inch  displacement  of  the  projector  directly  above  the 
viewer.  In  this  case,  the  projected  View  Window  subtends  angles  of  7CC  ver¬ 
tical  by  90°  horizontal  from  the  viewpoint.  On  this  View  Window,  there  is  an 
11  by  11  raster  pattern  which  forms  one  hundred  rectangular  blocks  as  depict¬ 
ed  on  Figure  A-l .  The  window  is  centered  on  the  surface  of  the  dome  at  the 
3-0  location  (240.0,  0.0,  0.0)  inches,  or  at  the  intersection  of  the-  X-axis 
with  the  dome  surface.  Table  A-l  shows  the  program  prompts  and  user  inputs 
to  describe  the  above  system.  Figures  A-2,  A-3,  and  A-4  show  the  required 
raster  shapes  to  be  projected  from  the  defined  projector  location  for  the  use 
of  projection  lenses  with  mapping  functions  of  Tan  e.  Sin  e,  and  e, 
respectively. 

The  second  projection/viewer  arrangement  is  a  projector  arrangement  with 
the  View  Window  and  projector  located  on  opposing  sides  of  the  dome  surface. 
In  this  system,  the  viewpoint  Is  again  coincident  with  the  center  of  a  20- 
foot  radius  dome.  However,  the  projector  is  located  behind  the  viewer  on  the 
surface  of  the  dome  at  the  3-D  location  (-218.0,  0.0,  101.0)  Inches,  and  the 
View  Window  Is  placed  at  the  3-D  location  (218.0,  0.0,  -101.0)  inches.  This 
arrangement  directs  the  optical  axis  of  the  projector  to  pass  through  the 
center  of  the  dome  on  Its  way  to  the  center  of  the  View  Window,  and  provides 
for  a  synmetrlcal  projection  onto  the  window  at  a  distance  of  twice  the  dome 
radius,  or  40  feet.  The  View  Window  contains  an  11  by  11  raster  pattern  and 
subtends  angles  of  160°  by  160°  from  the  viewpoint  as  shown  In  Figure  A-5. 
Figures  A-6,  A-7,  and  A-8,  with  Tables  A-6,  A-7,  and  A-8,  describe  the  re¬ 
quired  raster  shapes  for  Example  2.  Table  A-5  shows  the  necessary  Inputs  to 
describe  this  second  system. 
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Analysis  of  the  Sample  Systems 

Analysis  of  the  required  raster  shapes  for  a  projection  system  involves 
looking  at  the  graphic  mapping  outputs  and  their  corresponding  tabic  outputs. 
From  these  a  visualization  of  the  raster  shapes  is  obtained,  alcrg  with  the 
distortion  percentages  related  to  each  mapping. 

Examining  the  coordinate  tables  A-2,  A-3,  and  A-4,  which  are  for  the 
first  example,  specific  information  for  each  type  of  lens  ir.  this  projection 
arrangement  is  revealed.  The  tables  are  labeled  as  to  the  type  of  lens  map¬ 
ping  function  considered  along  with  identification  of  the  projector,  viewer, 
and  View  Window  locations.  The  tables  contain  the  coordinates  of  the  inter¬ 
section  points  for  the  distorted  raster  shape  and  the  coordinates  of  the 
non-disterted  Comparison  Plane.  Notice  that  the  coordinates  are  preceded  by 
the  labels  (X,  Y)  and  (XC,  YC)  which  denote  the  coordinates  of  the  view  plane 
and  Comparison  Plane,  in  turn.  These  coordinates  are  included  in  the  tables 

ir.  the  same  manner  that  the  raster  lines  are  drawn  on  the  View  Window;  that 

is,  left  to  right  and  top  to  bottom.  The  tables  provide  only  the  coordinates 
of  the  points  used  in  the  distortion  calculation  routine,  and  not  all  the 
points  of  the  11  by  11  pattern.  The  coordinates  of  all  the  points  are  avail¬ 
able,  but  are  not  included  here  due  to  the  size  of  the  tables. 

The  first  output  mapping  for  Example  1,  Figure  A-l,  is  a  tangent  mapping 
of  the  View  Window,  or  equivalently,  the  Comparison  Plane.  This  ic  the  raster 
shape  that  the  viewer  should  see  for  the  condition  of  no  distortion.  The 
tangent  mapping  for  this  projection  system.  Figure  A-2,  shows  a  maximum  dis¬ 
tortion  of  2.6  percent,  with  symmetrical  distortion  about  the  center  vertical 
line  of  the  raster  pattern.  All  the  distortion  percentages  for  this  mapping 
are  low  and  reflect  the  close  proximity  to  the  viewpoint  cf  a  projector  with  a 
F  Tan  e  mapping  function.  Looking  at  the  F  e  and  F  Sin  e  outputs,  the  maxi¬ 
mum  distortion  percentages  jump  to  15.7  percent  and  20.3  percent.  These 
percentages  are  large  and  reflect  the  departure  of  the  individual  lens  mapping 
functions  from  a  distortionless,  or  F  Tan  G,  lens.  In  this  system,  a  lens 
with  a  F  Tan  e  mapping  would  be  preferred  due  to  the  closeness  of  the  projec¬ 
tor  and  viewer,  and  the  less  than  90°  projection  angles  required. 
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The  second  example  is  offered  as  an  extreme  case  compared  to  the  preceding 
case.  The  separation  cf  the  projector  and  viewer  is  increased  to  2G  feet  and 
the  View  Window  is  required  to  fill  a  field  of  view  of  160°  by  160°  at  the 
viewpoint.  Examining  the  mappings  and  tables  for  this  second  system,  it  can 
be  seen  that  the  distorted  raster  shapes  for  all  the  mapping  functions  fall 
within  the  90°  reference  circle.  In  fact,  the  maximum  projection  angle  for 
any  of  the  lenses  is  approximately  80°.  This  angle  is  appropriate  considering 
the  View  Window  to  fill  a  160°  by  160°  field  of  view  at  the  viewpoint  with 
projection  from  a  distance  of  twice  the  dome  radius.  A  projection  lens  with 
a  total  field  of  view  of  90°  would  provide  more  than  enough  of  the  projection 
field  required  to  fill  the  View  Window. 

The  distortion  percentages  for  all  three  lens  types  are  quite  large  and 
are  approximately  equal.  Examining  Tables  A-6,  A-7,  and  A-8,  the  maximum 
distortion  required  for  the  F  e  mapping  is  slightly  greater  than  the  tangent 
mapping  (18.09  percent  vs.  18.90  percent),  while  the  F  Sin  e  mapping  requires 
the  most  raster  shaping  at  19.22  percent.  These  raster  shapes  may  be  hard  to 
implement  by  raster  shaping  alone  as  can  be  seen  in  Figures  A-6,  A-7,  and 
A-8. 

Other  information  to  be  gained  from  the  program  output  concerns  the  re¬ 
distribution  cf  the  raster  points  or.  the  projection  lens  target  plane.  In 
the  second  system,  the  points  are  crowded  together  near  the  edges  of  the 
target  plane  with  only  a  few  points  in  the  center  of  the  target  plane.  The 
effect  of  this  distribution  depends  on  the  projection  lens  in  use,  but  would 
tend  to  reduce  the  resolution  capacity  cf  the  imaging  system.  This  is  due  tc 
the  crowding  of  resolution  elements  {RESELS)  on  the  target  plane  where  most 
lenses  lose  resolution  capability,  and  lack  of  RESELS  where  most  lenses  have 
their  greatest  resolution. 
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IV  SUWARY 

The  program  MAPTAG  Is  a  very  useful  tool  for  determining  the  required  raster 
shape  to  be  projected  which  provides  for  distortionless  viewing.  The  graphic 
output  provides  visualization  of  distortions  encountered  In  spherical  screen 
projections.  The  table  output  allows  distortion  calculations  to  be  performed 
and  can  also  provide  for  equations  that  describe  the  raster  distortion,  line  by 
line.  The  subroutines  used  by  MAPTAG  have  been  written  In  general  form  to 
allow  their  use  in  building  other  specific  distortion  routines. 

Extension  of  the  program  to  provide  more  information  about  projection  sys¬ 
tems  is  possible  and  Is  being  Investigated.  These  programs  provide  an  excel¬ 
lent  basis  for  distortion  analysis  of  video  projection  systems  and  efforts 
are  being  taken  to  include  projection  of  generalized  imagery.  Additionally, 
the  type  of  screen  surface  used  in  the  program  is  being  expanded  to  include 
cylindrical  screens. 
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EXAMPLE  1 


RUN  MAPTAG 

DESCRIBE  PLANAR  INPUT  SCREEN 

ENTER  DESIRED  COORDINATES  FOR  PLACEMENT  OF  CENTER  OF  PLANE  IN  SPHERE 
COORDINATE  SYSTEM,  THE  X-AXIS  IS  POSITIVE  FORWARD  WITH  THE  Y-AXIS  POSITIVE 
TO  THE  LEFT,  AND  THE  Z-AXIS  IS  POSITIVE  UPWARDS 
240,0,0 

ENTER  VERTICAL  AND  HORIZONTAL  FOV'S  IN  DEGREES 
70,90 

ENTER  NUMBER  OF  POINTS  ALONG  VERTICAL  AND  HORIZONTAL  AXES 

11,11 

DESCRIBE  SPHERICAL  SCREEN  DISPLAY  SYSTEM 

ENTER  SPHERICAL  SCREEN  RADIUS 
240 

ENTER  PROJECTOR  POSITION  (FX.FY.FZ) 

0,0,12 

ENTER  VIEWPOINT  POSITION  (XE,YE,ZE) 

0,0,0 

ENTER  VERTICAL  AND  HORIZONTAL  OFFSET  FOR  PROJECTION  DIRECTION  (DEGREES) 

0,0 

ENTER  SCALE  FACTOR  ,  SCREEN  MAG  ,  YBIAS,  XBIAS 
NORMALLY  SCREEN  MAG  =  1 .  , XBIAS  =  0., YBIAS  =  0, 

TKE  SCALE  FACTOR  DETERMINES  THE  TFOV  (FSCALE  *  90) 

1,5, 1.1, 0.0 

PICK  WHICH  TYPE  OF  MAPPING  DESIRED,  FOR  A  TANGENT  MAPPING  ENTER  "1"  , 

FOR  AN  IDEAL  R-THETA  MAPPING  ENTER  "2"  ,  FOR  TANTHETA  MAPPING  WITH 
PRIMARY  DISTORTION  ,  ENTER  "3"  ,FOR  SIN  THETA  MAPPING,  ENTER  "4" 

1 

IF  TABLE  OUTPUT  IS  DESIRED  TYPE  TRUE  ;  IF  NOT  FALSE 
T 

TYPE  C  TO  CONTINUE  EXECUTION 

1 

$  C 

TO  CHANGE  VOFF, HOFF, MAG, FSCALE, LTABLE, TYPE  OF  MAPPING  ENTER  "T"  ,  IF  NOT 
ENTER  "F" 

F 

TO  CHANGE  SCREEN  MAGNIFICATION, ENTER  T  ,T0  STOP  ENTER  F 
F 

FORTRAN  STOP 

$ 


TABLE  A-l 
22 
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Figure  A1 
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TAN  THETA  MAPPING  COORDINATES 
THE  VIEWPOINT  IS  AT  : 


X 

=  0.000 

Y  -  0.000 

Z  =  0.000 

OBSERVING  A  LINEAR  RASTER  PATTERN  ON  THE  SCREEN  FILLING  A  FIELD  OF 
7G.000  DEGREES  VERTICALLY,  AND  90.000  DEGREES  HORIZONTALLY  ABOUT  A 

POINT  X= 

240.000  Y 

=  0.000  Z=  0.000  DEFINED  AS  CENTER  OF 

THE  PROJECTOR  IS  AT: 

X 

=  0. 

0000  Y  = 

O.OOC  Z  =  12.00C 

THE  PROJECTION  DIRECTION  IS  OFFSET  BY  V=  0.000  H  =  0.000  FROM  ' 

DEFINED 

FOV  CENTER. 

THE 

RADIUS  OF  ThE  SPHERICAL  SCREEN  IS  240. COO 

X 

1.033 

0.000 

-1.033 

Y 

0.693 

0.712 

0.693 

XC 

1.000 

0.000 

-1.000 

YC 

0.700 

0.700 

0.700 

X 

0.998 

0.000 

-0.998 

Y 

-0.021 

0.000 

-0.021 

XC 

1.000 

0.000 

-1.000 

YC 

O.OOC 

0.000 

0.000 

X 

G.S64 

C.000 

0.964 

Y 

-0.702 

-0.685 

-0.702 

XC 

1.000 

0.000 

-1 .000 

YC 

-0.700 

-0.700 

-0.700 

DISTORTION  PERCENTAGES 

THE  TARGET  PLANE  HEIGHT 

IS  1.397 

1 

2 

3 

2.59$ 

0.96$ 

2.59$ 

4 

5 

6 

1.48$ 

0.00$ 

1 .48$ 

7 

8 

9 

2.44$ 

0.96$ 

2.44$ 

TABLE  A-2 
24 
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R-THETA  MAPPING  COORDINATES 
THE  VIEWPOINT  IS  AS  : 

X  =  0.0000  Y  =  0.000  Z  «  0.000 

OBSERVING  A  LINEAR  RASTER  PATTERN  ON  THE  SCREEN  FILLING  A  FIELD  OF 
70.000  DEGREES  VERTICALLY,  AND  90.000  DEGREES  HORIZONTALLY  ABOUT  A 
POINT  X=  240.000  Y=  0.000  Z=  0.000  DEFINED  AS  CENTER  OF  "FOV". 

THE  PROJECTOR  IS  AT  : 

X  =  0.0000  Y  =  0.000  Z  =  12.000 

THE  PROJECTION  DIRECTION  IS  OFFSET  BY  V=  0.000  R  =  0.000  FROM 
THE  DEFINED  FOV  CENTER. 

THE  RADIUS  OF  THE  SPHERICAL  SCREEN  IS  240.000 


X 

0.742 

O.OOC 

-0.742 

Y 

0.498 

0.619 

0.498 

XC 

1.000 

0.000 

-1.000 

YC 

0.700 

0.700 

0.700 

X 

0.784 

C.000 

-0.784 

Y 

-0.016 

O.OOC 

-0.016 

XC 

l.COO 

0.000 

-1.000 

YC 

O.OOC 

0.000 

O.OGO 

X 

0.706 

C.000 

-0.706 

Y 

-0.514 

-0.601 

-0.514 

XC 

1.000 

0.000 

-1.000 

YC 

-0.700 

-0.700 

-0.700 

DISTORTION 

PERCENTAGES 

THE 

TARGET  PLANE  HEIGHT 

IS  l.; 

1 

2 

3 

13.97* 

0.74% 

13.97% 

4 

5 

6 

7.22% 

0.00% 

7.22% 

7 

8 

9 

15.66% 

0.74% 

15.66% 

TABLE  A-3 
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SIN  THETA  MAPPING  COORDINATES 


THE  VIEWPOINT  IS  AT  : 

X  ■  O.OCCO  Y  =  0.000  Z  =  0.000 

OBSERVING  A  LINEAR  RASTER  PATTERN  ON  THE  SCREEN  FILLING  A  FIELD  OF 
70.000  DEGREES  VERTICALLY,  AND  9C.000  DEGREES  HORIZONTALLY  ABOUT  A 
POINT  X=  240. COO  Y=  0.000  Z=  0.000  DEFINED  AS  CENTER  OF  "F0V\ 

THE  PROJECTOR  IS  AT  : 

X  *  0.0000  Y  =  0.000  Z  =  12.000 

THE  PROJECTION  DIRECTION  IS  OFFSET  BY  V=  C.000  H  =  0.000  FROM 

THE  DEFINED  FOV  CENTER. 


THE  RADIUS  OF  THE  SPHERICAL  SCREEN  IS  240. COO 


X 

0.647 

0.000 

-0.647 

Y 

0.434 

0.580 

0.434 

XC 

1.000 

0.000 

-1.000 

YC 

0.700 

0.700 

0.700 

X 

0.706 

0.000 

-0.706 

Y 

-0.015 

C.000 

-0.015 

XC 

1.000 

0.000 

-1.000 

YC 

0.000 

0.000 

0.000 

X 

0.619 

G.GOC 

-0.619 

Y 

-0.451 

-0.565 

-0.451 

XC 

1.000 

0.000 

-1 .000 

YC 

-0.700 

-0.700 

-0.700 

DISTORTION 

PERCENTAGES 

THE 

TARGET  PLANE  HEIGHT 

IS 

1 

2 

3 

19.17% 

0.64% 

19.17% 

4 

5 

6 

9.82% 

0.00% 

9.82% 

7 

8 

9 

20.32% 

0.64% 

20.32% 

TABLE  A-4 
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EXAMPLE  2 

SAMPLE  RUN  OF  HAPTAG 

RUN  MAPTAG 

DESCRIBE  PLANAR  INPUT  SCREEN 

ENTER  DESIRED  COORDINATES  FOR  PLACEMENT  OF  CENTER  OF  PLANE  IN  SPHERE 
COORDINATE  SYSTEM  .THE  X-AXIS  IS  POSITIVE  FORWARD  WITH  THE  Y-AXIS  POSITIVE 
TO  THE  LEFT,  AND  THE  Z-AXIS  IS  POSITIVE  UPWARDS 

216,0,-101 

ENTER  VERICAL  AND  HORIZONTAL  FOVS  IN  DEGREES 
160,160 

ENTER  HUMBER  OF  POINTS  ALONG  VERTICAL  AND  HORIZONTAL  AXES 

11,11 

DESCRIBE  SPHERICAL  SCREEN  DISPLAY  SYSTEM 
ENTER  SPHERICAL  SCREEN  RADIUS 


240 

ENTER  PROJECTOR  POSITION  (FX.FY.FZ) 
-216,0,101 

ENTER  VIEWPOINT  POSITION  (XE.YE.ZE) 


0,0,0 

ENTER  VERTICAL  AND  HORIZONTAL  OFFSET  FOR  PROJECTION  DIRECTION  (DEGREES) 


0,0 

ENTER  SCALE  FACTOR  ,  SCREEN  MAG  .  YBIAS,  XBIAS  NORMALLY  SCPEEN  MAG  =  1. 
.XEIAS  *  C.  , YEIAS  =  0.  THE  SCALE  FACTOR  DETERMINES  THE  TFOV  (FSCALE 
*  90) 

1.5, 1,0,0 

PICK  WHICH  TYPE  OF  MAPPING  DESIRED,  FOR  A  TANGENT  MAPPING  ENTER  "1"  ,  FOR 
AN  IDEAL  R-THETA  MAPPING  ENTER  "2"  ,  FOR  TANTHETA  MAPPING  WITH  PRIMARY 
DISTORTION  ,  ENTER  "3"  .FOR  SIN  THETA  MAPPING,  ENTER  "4" 


IF  TABLE  OUTPUT  IS  DESIRED  TYPE  TRUE  ;  IF  NOT  FALSE 
T 

TYPE  C  TO  CONTINUE  EXECUTION 

1 

$  C 

TO  CHANGE  VOFF, HOFF, MAG, FSCALE .LTABLE .TYPE  OF  MAPPING  ENTER  "T"  ,  IF  NOT 
ENTER  "F" 

T 

ENTER  VERTICAL  AND  HORIZONTAL  OFFSET  FOR  PROJECTION  DIRECTION  (DEGREES) 
O.C 

ENTER  SCALE  FACTOR  ,  SCREEN  MAG  ,  YBIAS,  XBIAS 
NORMALLY  SCREEN  MAG  =  1 .  .XBIAS  *  0. .YBIAS  -  C.  THE  SCALE  FACTOR 
DETERMINES  THE  TFOV  (FSCALE  *  90) 
l.S, 1,0,0 

PICK  WHICH  TYPE  OF  MAPPING  DESIREO,  FOR  A  TANGENT  MAPPING  ENTER  "1"  , 

FOR  AN  IDEAL  R-THETA  MAPPING  ENTER  "2"  ,  FOR  TANTHETA  MAPPING  WITH 
PRIMARY  DISTORTION  ,  ENTER  ”3"  ,FOR  SIN  THETA  MAPPING,  ENTER  "4" 


IF  TABLE  OUTPUT  IS  DESIRED  TYPE  TRUE  ;  IF  NOT  FALSE 
F 

NO  OUTPUT  FOR  TABLE 
TYPE  C  TO  CONTINUE  EXECUTION 

1 

$  C 

TO  CHANGE  VOFF.HOFF.MAG, FSCALE, LTABLE, TYPE  OF  MAPPING  ENTER  "T"  ,  IF  NOT 
ENTER  "F" 

F 

TO  CHANGE  SCREEN  MAGNIFICATION, ENTER  T  ,T0  STOP  ENTER  F 
F 

FORTRAN  STOP 

i 


TABLE  A-5 
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COMPARISON  PLANE 
VIEW  WINDOW 
160°  x  160° 


IH  Report  332 


NAVTRAEQUIPCEN  IH-332 


TAN  THETA  MAPPING  COORDINATES 


THE  VIEWPOINT  IS  AT  : 

X  =  0.0000  Y  *  0.000  Z  =  0.000 

OBSERVING  A  LINEAR  RASTER  PATTERN  ON  THE  SCREEN 
FILLING  A  FIELD  OF  160.000  DEGREES  VERTICALLY, 

AND  160.000  DEGREES  HORIZONTALLY 
ABOUT  A  POINT  X=  218.000  Y=  0.000 
Z=  -101.000  DEFINED  AS  CENTER  OF  "FOV". 

THE  PROJECTOR  IS  AT  : 

X  =  -218.0000  Y  =  0.000  Z  =  101.000 

THE  PROJECTION  DIRECTION  IS  OFFSET  BY  V=  0.000  H 
FROM  THE  DEFINED  FOV  CENTER. 


THE  RADIUS  OF  THE  SPHERICAL  SCREEN  IS 


X 

0.624 

0.000 

-0.624 

Y 

0.624 

0.838 

0.624 

XC 

5.671 

0.000 

-5.671 

YC 

5.671 

5.671 

5.671 

X 

0.838 

0.000 

-0.838 

Y 

0.000 

0.000 

0.000 

XC 

5.671 

0.000 

-5.671 

YC 

0.000 

0.000 

0.000 

X 

0.624 

0.000 

-0.624 

Y 

-0.624 

-0.838 

-0.624 

XC 

5.671 

0.000 

-5.671 

YC 

-5.671 

-5.671 

-5.671 

DISTORTION 

PERCENTAGES 

1 

2 

3 

18.092 

0.002 

18.092 

4 

5 

6 

0.0C2 

0.002 

0.002 

240.000 


7  8  9 

18.092  0.002  18.092 


TABLE  A-6 
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R-THETA  MAPPING  COORDINATES 
THE  VIEWPOINT  IS  AT  : 

X  =  0.0000  Y  =  0.000  Z  =  0.000 

OBSERVING  A  LINEAR  RASTER  PATTERN  ON  THE  SCREEN 
FILLING  A  FIELD  OF  160.000  DEGREES  VERTICALLY, 

AND  160,000  DEGREES  HORIZONTALLY 
ABOUT  A  POINT  X=  218.000  Y=  0.000 

Z=  -101.000  DEFINED  AS  CENTER  OF  "FOV". 

THE  PROJECTOR  IS  AT  : 

X  =  -218.000  Y-  C.000  Z  =  101. OCO 

THE  PROJECTION  DIRECTION  IS  OFFSET  BY  V=  0.000  H  = 
FROM  THE  DEFINED  FOV  CENTER. 

THE  RADIUS  OF  THE  SPHERICAL  SCREEN  IS  240.000 


X 

0.511 

o.oco 

-0.511 

Y 

G.511 

0.698 

0.511 

XC 

5.671 

0.000 

-5.671 

YC 

5.67. 

5.671 

5.671 

X 

0.C98 

0.000 

-0.698 

Y 

0.000 

0.000 

0.000 

XC 

5.671 

0.000 

-5.671 

YC 

C.COO 

0.000 

0.000 

X 

0.511 

0.000 

-0.511 

Y 

-0.511 

-0.698 

-0.511 

XC 

5.671 

0.000 

-5.671 

YC 

-5.671 

-5.671 

-5.671 

DISTORTION  PERCENTAGES 


1  2  3 


18.902 

0.00°/ 

18.902 

4 

5 

6 

0.002 

r'.C01 

0.CO2 

7 

S 

9 

IS.  SO0/ 

0.002 

18.902 

TABLE  A-7 
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SIN  THETA  MAPPING  COORDINATES 
THE  VIEWPOINT  IS  AT  : 

X  =  0.0000  Y  =  0.000  7  =  0.000 

OBSERVING  A  LINEAR  RASTER  PATTERN  ON  THE  SCREEN 
FILLING  A  FIELD  OF  160.000  DEGREES  VERTICALLY, 

AND  160,000  DEGREES  HORIZONTALLY 
ABOUT  A  POINT  X=  218.000  Y=  0.000 

Z=  -101. OOC  DEFINED  AS  CENTER  OF  "FOV". 

THE  PROJECTOR  IS  AT  : 

X  =  -218.000  Y=  0.000  Z  =  101.000 

THE  PROJECTION  DIRECTION  IS  OFFSET  BY  V=  0.000  H  = 
FROM  THE  DEFINED  FOV  CENTER. 

THE  RADIUS  OF  THE  SPHERICAL  SCREEN  IS  240. COO 


X 

0.468 

0.000 

-0.468 

Y 

0.468 

0.642 

0.468 

XC 

5.671 

0.000 

-5.671 

YC 

5.671 

5.671 

5.671 

X 

0.642 

0.000 

-0.642 

Y 

0.000 

0.000 

0.000 

XC 

5.671 

0.000 

-5.671 

YC 

0.000 

0.000 

0.000 

X 

0.468 

C.000 

-0.468 

Y 

-0.468 

-0.642 

-0.468 

XC 

5.671 

0.000 

-5.671 

YC 

-5.671 

-5.671 

-5.671 

DISTORTION  PERCENTAGES 

1 

2 

3 

19.22% 

0.00% 

19.22% 

4 

5 

6 

0.00% 

0.00% 

0.00% 

7 

8 

9 

19.22% 

0.00% 

19.22% 

TABLE  A -8 
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110°  FOV 


SIN  THETA  MAPPING 


< 


NAVTRAEQU I PCEN  IH-332 

CONVERSION  FRON  F  TAN  e  KAPPINGS  TO  Fg,  F  Sin  6, 
AND  F  TAN  G  WITH  PRIMARY  DISTORTION  MAPPINGS 


F  Tan  e  mapping 

Using  spherical  angles  4-  and  e  tc  define  angles  along  the  X  and  Y  axes  of 
the  mapping  plane  relative  tc  the  origin  (center  of  plane)  for  a  point  cn  the 
plane,  the  coordinates  of  the  point  are  found  by  the  equations: 


Xt  =  F  *  Tan  $ 
Yt  =  F  *  Tan  0 

F  e  mapping 


A  point  on  the  F  Tan  £  napping  plane  (for  F  =  1)  is  located  at  a  radial 
distance  from  the  origin  given  by: 

Rt  =  [(Tan  *)2  +  (Tan  0)2]l/z. 

Considering  the  difference  between  the  radial  distance  for  a  Tan  o  mapping 
and  a  g  mapping  for  the  same  point  in  the  field,  the  e  mapping  radial  distance 
is  given  by: 


Re  = 


tan  a’ 


or  equivalently. 


Where  CR  Is  the  conversion  ratio  constant  for  a  particular  point  cn  the  plane 
(Tan  e),  a  Is  the  angle  from  the  optical  axis. 

To  find  the  X  and  Y  coordinates  on  the  e  mapping  plane,  expand  the  conver¬ 
sion  equation  to: 

C(xe)j  +  (Ye)23l/2  *  [(Xt)2  +  <rt)2]‘/2*  CR  . 
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and  using  the  property  of  a  scaler  multiplication  of  a  vector  the  resulting 
coordinates  for  a  F  0  mapping  are  given  by: 


xe  -  xt  *  CR  . 


Ye  "  Yt  *  CR  * 


F  sin  e  Mapping 


Again  considering  the  radial  distance  for  a  point  on  a  Tan  e  mapping, 
if  the  radial  distance  (Rs)  for  a  Sin  o  mapping  is  desired,  then 


R 


s 


*  Sin  a 

Tan 'a  * 


or. 


[(Xs)2  ♦  (YSP]1/2  •  «xt)2  ♦  (Yt)2]>/2  X  4aS-i-  • 


using  the  scaler  multiplication  process  again. 


X 


s 


y  £  S 1  0  Cl 

- 

Tan  a  , 


Y 


s 


_  w  x  S in  ct 
Yt  Tan  a  . 


Tan  0  mapping  with  Primary  Distortion 

Primary  distortion  of  a  Tan  e  lens  is  a  distortion  of  the  radial  position 
of  a  point  on  the  mapping  plane  relative  to  the  square  of  the  radial  distance 
to  that  point  multiplied  by  a  constant  factor  (DF),  i.e.: 

Rp  •  Rt  +  (DF  *  Rt2)  , 


or. 
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Rp  “  Rt  *  (  1  +  (OF  *  Rt)). 

Using  the  general  formula  for  radial  distance  the  above  equation  expands 
to: 

[(Xp)2  +  (Yp)2]i/2  =  [(Xt)2  +  (Yt)2]i/2  Ml  +  DF  *  Rt)), 
again  scalar  multiplication  process  leads  to: 
xp  =  Xt  *  (  1  +  (DF  *  RtJ)  , 

Yp  •  Yt  M  1  ♦  (CF  +  Rt))  , 

where  Xp  and  Yp  are  the  coordinates  of  the  distorted  mapping.  Obviously, 
letting  CF  =  0  will  give  a  Tan  o  mapping,  with  a  negative  DF  value  giving 
barrel  distortion  and  a  positive  DF  value  giving  pincushion  distortion. 

The  figures  below  illustrate  the  difference  between  the  mappings  through 
a  lens  of  a  rectangular  object  for  the  four  general  mappings  described  above. 
Notice  that  for  the  proper  choice  of  DF,  an  F  c  mapping  or  F  Sin  e  mapping 
can  be  approximated  if  the  field  angles  are  r.ot  too  large. 
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Y 


Figure  B-l.  Relative  Distortions  of  a  Rectangular  Object 
by  Various  Lens  Happing  Functions. 
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PROGRAM  FLOWCHARTS  AND  CODING 
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MAPTAG 
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|  STORE  CORNER  POINT  (COMPARISON  POINT) 

_ i _ 


CALCULATE  ROTATION  CONSTANTS  : 
TO  ROTATE  INPUT  PLANE  TO  BE 
PERPENDICULAR  TO  LINE  OF  SIGHT 
TO  ROTATE  POINTS  ON  SPHERE  TO 
CENTER  PCTR  ON  OPTICAL  AXIS 


-  LOOP  FOR  2nd  TO  NPT  INPUT  POINTS 


CALL 

SPTFANG 

STORE  POINTS 

ON  SPHERE 

TILT,  ROTATE,  OFFSET  INPUT  PLANE 
TO  FORM  VIEW  WINDOW 


FIND  SPHERICAL  ANGLES  FROM  VIEWPOINT 
TO  POINTS  ON  VIEW  WINDOW 


PROJECT  VECTORS  THROUGH  VIEW  WINDOW 
ONTO  SPHERE  FROM  VIEWPOINT 


FIND  PERSPECTIVE  VIEW  OF  POINTS  ON 
SPHERE  FROM  PROJECTION  POINT 
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points 


CLEAR 

GRAPHIC 

SCREEN 
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INPUT: 


OUTPUT: 


SUBROOTINE  TROPLN 


t 

vi. 


RETURN 


J 


ROTATE  PLANE  ABOUT  Y-AXIS 
BY  KV  DEGREES. 


ROTATE  PLANE  ABOUT  Z-AXIS 
BY  RH  DEGREES. 


POINT  CM  PLANE/P (PX,PY,PZ)  ,  VERTICAL  ROTATION  ANGLE  (KV)  ,  HORIZONTAL 
ROTATION  ANGLE  (RH)  ,  PLANE  CENTER  PCTR(PXC,PYC,PZC)  ,  ORIGIN  OF  PLANE 
(XO,YO,ZO) . 


POINT  ON  PLANE  P(PX,PY,PZ)  -  (TILTED, ROTATED  AND  OFFSET) 
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corn 

00C2 


0003 

0004 

0005 


C006 


0007 


oc  : 

0009 

C010 


0011 

0012 


22-JAN-1S81  15:45:10  VAX  11  FORTRAN  IV-FLUS  VI. 2-13  Paae  1 

SUBDIST. FOR. 1 

C 

C 

£  ★**★★★**★* 

C  *  * 

C  *  TROPLN  * 

C  *  * 

£  **★★★★★★*■* 

C 

C  THIS  SUBROUTINE  COVERTS  POINTS  ON  A  Y-Z  PLANE  CENTERED  AT 
C  AN  ORIGIN  TO  POINTS  ON  A  PLANE  THAT  IS  TILTED  ABOUT  THE  :-A ■ . 

C  EY  RV  DEGREES  AND  THEN  ROTATED  Rh  DEGREES  AEC'UT  THE  Z-/>  ’ ' 

C  THE  PLANE  IS  THEN  TRANSLATED  FROM  THE  ORIGIN  TO  ( F > C , F  • 0 , :  . 

C  THE  PLANE  PRODUCED  WILL  EE  TANGENT  TC  THE  PROJECTION 

C  DIRECTION  IF  CTHE  ROTATION  ANGLES  ARE  THE  NEGATIVE  OF  THE 

C  ANGLES  FROM  PROJECTOR  TO  (PXC.PYC.PZC)  OF  THE  SPHERE  CENTtF r: 
C  COORDINATE  SYSTEM. 

C 

C 

SUBROUTINE  TRCPLN  (RV,RH,PXC,PYC,PZC,XC,Y0,70,PX,PY,P7) 
IMPLICIT  DOUBLE  PRECISION  (D,F,P,P,X,Y,Z) 

C 

C 

C  TRANSLATE  PLANE  TC  ORIGIN  AT  (XO.YO.ZO) 

C 

PX  =  PX  -  XO 
PY  =  PY  -  YC 
PZ  =  PZ  -  ZO 
C 
C 

C  ROTATE  PLANE  ABOUT  Y-AXIS  DY  RV  DEGREES 
C 

CALL  VROTAT(RV,PX,PY,PZ) 

C 

C  ROTATE  ABOUT  Z-AXIS  BY  RH  DEGREES 
C 

CALL  HROTAT(RH ,PX»PY,PZ) 

r 

o 

C 

C  TRANSLATE  PLANE  FROM  PROJECTION  POINT  TC  POSITION  CENTERED  AT 
C  PXC,PYC,PZC. 

C 

PX  «  PX  +  PXC  -  XC 
PY  =  PY  +  PYC  -  YC 
PZ  =  PZ  +  PZC  -  ZC 
C 
C 

RETURN 

END 
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INPUT  : 
OUTPUT: 


SUBROUTINE  SANGFPT 


START 


+/-  90  degrees 


+/-  180  degrees 


+/~  180  degrees 


DEGREES  TO  RADIAN  FACTOR  (PI)  ,  POINT  F(FX,FY,FZ)  ,  POINT  P(PX,PY,PZ)  . 

'VANG'  (VERTICAL  +/“  90  DEGREES  FROM  HORIZON)  ,  'TILT'  (  -'VANG'),  'VANGR' 
(VERTICAL  +/-  180  DEGREES) ,  'HANGR'  (  HORIZONTAL  +/-  180  DEGREES) , 

'DIST'  (  DISTANCE  FF  ) 


5C 
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30- JAN-1 981  12:12:16  VAX-11  FORTRAN  IV-PLUS  VI. 2-13  Page  7 

TESTSUB.FOR.l 


C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


************ 

*  * 

*  SANGFPT  * 

*  * 

************ 

THIS  SUBROUTINE  FINDS  SPHERICAL  ANGLES  TO  POINT  PX.PY.PZ 
FROM  POINT  FX.FY.FZ.  ANGLES  ARE  EXPRESSED  IN  RADIANS.  ALSO 
THE  DISTANCE  FROM  F  TO  P  IS  RETURNED  IN  THE  SAME  UNITS  AS 
POINT  LOCATIONS  ARE  GIVEN. 

THREE  ANGULAR  SIGN  CONVENTIONS  ARE  USED  IN  THIS  PROGRAM 

1.  POSITIVE  ABOVE  HORIZON  ,  NEGATIVE  BELOW  HORIZON 
CONVENTION  .  VERTICAL  ANGLES  ONLY  (+/-  9C  DEG.) 

VALUES  RETURNED  AS  "VANG". 

2.  POSITIVE  FOR  COUNTER-CLOCKWISE  ROTATION  ,  NEGATIVE 
FOR  CLOCKWISE  ROTATION  AS  VIEWED  FROK  POSITION  ON 
A  POSITIVE  AXIS  LOOKING  TOWARDS  PLANE  DEFINED  BY 
REMAINING  AXES.  VERTICAL  AND  HORIZONTAL  ANGLES; 

VALUES  RETURNED  AS  1 VANGR'  AND  'HANGR' 

3.  VALUE  FOR  TILTING  OF  INPUT  PLANE  BY  TROPLN  SUBROUTINE 
THIS  ANGLE  IS  NEGATIVE  OF  VANG. 


001 

r 

SUBROUTINE  SANGFPT(PI ,PX,PY ,PZ ,FX ,FY ,FZ 
TILT  1) 

C02 

L 

r 

IMPLICIT  DOUBLE  PRECISION  (A,B,C,D,F,H, 

L 

c 

FIND  DISTANCE  FROM  POINT  F  TO  POINT  P 

003 

A  =  PX  -  FX 

004 

B  =  PY  -  FY 

OOS 

C  =  PZ  -  FZ 

ooe 

r 

DIST  *  DSQRT(A  *A  +  E*B  +  C*C) 

c 

r 

FIND  DIRECTION  COSINES  TO  P  FROM  F 

007 

L 

DYC  =  B  /  DIST 

008 

r 

DZC  =  C  /  DIST 

L 

c 

r 

FIND  SPHERICAL  ANGLES  TO  P  FROM  F 

009 

L 

VANG  =  DAS  IN (DZC) 

c 

STORE  ABSOLUTE  VERTICAL  ANGLE  (+/-  180 

010 

IF  (A  .LE.  0.  .AND.  C  .GT.  0. )THEN 

on 

VANGR  *  -  18C.*PI  +  VANG 

012 

ELSE 

013 

IF  (  .LT.  0.  .AND.  C  .LT.  0.)  THEN 

014 

VANGR  =  18C.*PI  +  VANG 

01E 

ELSE 

016 

VANGR  ■=  -  VANG 

017 

ENDIF 

018 

ENDIF 
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TESTSUB.FOR.1 

C  STORE  TITLE  VALUE  (+/-  90  DEG.) 

C 

0019  TILT  *  -  VANG 

C 
C 

C  FIND  PROJECTION  OF  VECTOR  ONTO  X-Y  PLANE 
C 


0020 

r 

CS  =  DSQRT { 1 .0  -  DZC*DZC) 

0021 

w 

IF  (CS  .EQ.  0.0)  THEN 

0022 

HANGR  *  0.0 

0023 

ELSE 

0024 

ARGH  =  DYC  /  CS 

0025 

IF  (  ARGH  .GT.  1.0)  ARGH  =  1.0 

0026 

D 

IF  (  ARGH  .LT.-1.0)  ARGH  *-1.0 

WRITE(2,*) , 'ARGH  ='  CS  *',CS 

0027 

HANGR  =  DASIN(ARGH) 

r 

ENDIF 

L 

c 

CONVERT  HORIZONTAL  ANGLES  TO  +/-  180  DEGREE  CONVENTION 

0029 

l 

IF  (  B  .LE.  0.  .AND.  A  ,LT.  0.  )  HANGR  *  -180.  *  PI  -  HANGR 

0030 

r 

IF  (  B  .GT.  0.  .AND.  A  .LT.  0.  )  HANGR  *  180.  *  PI  -  HANGR 

0031 

RETURN 

0032 

END 

NAVTRAEQU I PC EN  IH-332 


SUBROUTINE  SPTFANG 


INPUT:  RADIUS  (R)  ,  VERTICAL  ANGLE  (VANG)  ,  HORIZONTAL  ANGLE  (HANG) , 
PROJECTION  POINT  F(FX,FY,FZ) . 

OUTPUT:  INTERSECTION  POINT  ON  SCREEN  T(TX,TY,TZ)  . 
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TESTSUB.FOR.l 

C 

C 

C  *********** 

c  *  * 

C  *  SPTFANG  * 

C  *  * 

Q  *********** 

C 

c 

C  THIS  SUBROUTINE  FINDS  THE  POINT  OF  INTERSECTION  OF  A  VECTOR 

C  PROJECTED  FROM  A  POINT  (  FX,FY,FZ  )  AT  SPHERICAL  ANGLES  OF  'VANG’ 

C  AND  'HANG'  ONTO  THE  SURFACE  OF  A  SPHERE  WITH  THE  CENTER  DEFINED 

C  AS  THE  ORIGIN  OF  THE  3-D  SYSTEM.  THE  COORDINATE  SYSTEM  OF  THE 

C  PROJECTION  POINT  IS  PARALLEL  TO  THE  SCREEN  COORDINATE  SYSTEM. 

C  THE  POINT  RETURNED  IS  GIVEN  BY  (  PX,PY,PZ  )  AND  IS  ON  THE 
C  SURFACE  OF  THE  DESCRIBED  SPHERE  OF  RADIUS  (R) 

C 

OCOl  SUBROUTINE  SPTFANG  ( R, VANG, HANG, FX, FY ,FZ,TX,TY,TZ) 

0002  IMPLICIT  DOUBLE  PRECISION  (D.F.H.R.T.V) 

C 

C  FIND  DIRECTION  COSINES  FROM  F  TO  T 

C 

0003  DCV  =  DCOS(VANG) 

0004  DZP  =  DSIN(VANG) 

OOC5  DYP  =  DCV  *  DSIN(HANG) 

0006  DXP  =  DCV  *  DCOS(HANG) 

C 

C  FIND  LENGTH  OF  DISPLACEMENT  VECTOR  DUE  TO  PROJECTOR  OFF  CENTER 
C  OF  SPHERE 

C 


0007 

r 

DVMS  =  FX  *  DXP  +  FY  *  DYP  +  FZ  *  DZP 

C 

c 

r 

FIND  DISTANCE  FROM  PROJECTION  POINT  TO  POINT  ON  SCREEN 

0008 

l 

r 

DVAL  =  DABS  (DVMS*DVKS  +  R*R  -  FX*FY  -FY*FY  -  FZ*FZ) 

0009 

l 

r 

DIST  =  -  DVMS  +  DSQRT (  DVAL  ) 

t- 

c 

r 

FIND  COORDINATES  OF  POINT 

0010 

L 

TX  =  FX  +  DIST  *  DXP 

0C11 

TY  =  FY  +  DIST  *  DYP 

0012 

TZ  =  FZ  +  DIST  *  DZP 

0013 

t 

RETURN 

0014 

END 
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INPUT 


OFFSET 


SUBROUTINE  PSPECTIVE 


VIEWPOINT  F(FX,FY,FZi  ,  POINT  TO  BE  VIEWED  T(TX,TY,TZ) ,  ANGLES  TO 
ALICN  \  AXIS  TO  PCV  CENTER  (  VOFF  ,H3FF)  ,  ANGLES  TO  OFFSET  X  AXIS 
(  VCFFP,HOFFP) ,  MAPTYPE,  PRIMARY  DISTORTION  FACTOR  (  DFACTOR  ) . 

COORDINATES  OF  POINT  ON  VIEWPLANE  (  XOUT.YOUT  ) . 
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0001 


00C2 


0003 

0004 

0005 


0006 

0007 


0008 

0009 


C 

C 

c 

c 

c 

c 

c  *  * 

C  *  PERSPECTIVE  * 

C  *  * 

C  ************************** 

c 

C  THIS  SUBROUTINE  IS  USED  TO  FIND  THE  PERSPECTIVE  MAPPING  OF 
C  POINTS  IN  3-D  SPACE  WHERE  THE  COORDINATE  SYSTEM  IS  DESCRIBED 

C  AS:  X  IS  POSITIVE  FORWARD,  Y  IS  POSITIVE  TO  THE  RIGHT,  AND 

C  Z  IS  POSITIVE  UPWARDS  (  RIGHT-HAND  SYSTEM  ). 

C  THE  PERSPECTIVE  IS  FOUND  FROM  A  POINT  (FX.FY.FZ)  WITH  THE 
C  VIEWING  VECTOR  OFFSET  FROM  OH.OV  DEGREES  BY  -VOFF  DEGREES 
C  AND  -HOFF  DEGREES  IN  AN  ANGULAR  SYSTEM  THAT  DEFINES  THE 
C  DIRECTION  OF  PCSTIVE  ROTATION  AS  COUNTER-CLOCKWISE  ROTATION 
C  OF  TWO  AXES  WHEN  VIEWED  FROM  THE  REMAINING  AXIS. 

C 

C  Z 

C  ' 

C  * 

C  ’ 

C  ' 


C  OV.OH  X< - *»  EYE 

C  / 

C  / 

C  / 

C  v' 

C  V  0V.90H 

C  Y 

C 


C  THE  POINT  TO  BE  MAPPED  IS  (TX.TY.TZ)  ACCORDING  TO  MAPTYPE 
C  (1-TAN-THETA,2-THETA,3-TAN-THETA  W /  DFACTOR, 4-SIN-THETA) 

C 

C  VOFFP  AND  HOFFP  ALLOW  FURTHER  OFFSETTING  OF  THE  VIEWING  AXIS 

C  XOUT  AND  YOUT  ARE  THE  COORDINATES  FOR  PLANAR  PERPECTIVE  VIEW 

C 
C 

SUBROUTINE  PSPECTIVE(FX,FY,FZ,TX,TY,TZ,VOFF,HOFF,VOFFP,HOFFP 
1 MAPTYPE .DFACTOR , XOUT , YOUT) 

C 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

C 

C  TRANSLATE  ORIGIN  TO  VIEWPOINT  POSITION  (FX.FY.FZ) 

TX  =  TX  -  FX 
TY  =  TY  -  FY 
TZ  =  TZ  -  FZ 
C 

C  POINT  Xv  AXIS  TOWARD  FOV  CENTER  AND  FIND  NEW  COORDINATES 
C 

CALL  H ROT AT  (HOFF,TX,TY,TZ) 

CALL  VROTAT  (VOFF, TX.TY.TZ) 

C 

C  OFFSET  VIEWING  AXIS  IF  HOFFP  AND  VOFFP  ARE  NOT  ZERO 
C 

IF( VOFFP  .NE.  0.0)  CALL  VROTAT  (VOFFP, TX.TY.TZ) 

IF(HOFFP  .NE.  0.0)  CALL  HROTAT  (HOFFP .TX.TY.TZ) 


C 
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0010 

0011 

0012 

0013 

0014- 


0015 

0016 


0017 

0018 

0019 

0020 

0021 


0022 

0023 

0024 

0025 

0026 

0027 


0028 

0029 

0030 

0031 

0032 

0033 
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C  FIND  TANGENTS  OF  ANGLES  TO  POINTS 

C 

C 

IF  (TX  .EQ.  0.)  THEN 
C  SET 

XOUT  =  0. 

TOOT  -  0. 

RETURN 

ENDIF 

C 

C 

C 

TTHETA  -  TZ  /  TX 
TANPHI  -  TY  /  TX 
C 
C 
C 

CMfmmmtnnmftmmmHwmmnmmmtffmmfmmti 

ciminiiniimimniminnimiiimimimniniinmniinni 

TF{  MAPTYPE  .EQ.  1)  THEN 
C  FIND  TANGENT  PROJECTIONS  ONTO  PLANE 
C 

XOUT  =  TANPHI 
YOUT  »  TTHETA 
C 

RETURN 

C 

ENDIF 

ciinininimnimiinniiiiniimniiinnniimimnmnmni 

c 

c 

C  ELSE  FIND  TANGENT  RADIAL  DISTANCE  (TRDIST)  AND  CHECK  FOR  MAPTYPE 
C 

TRDIST  =  DSQRT( TANPHI  *  TANPHI  +  TTHETA  *  TTHETA) 

C 

IF(TRDI5T  .EQ.  0.)  THEN 
XOUT  =  0. 

YOUT  =  0. 

C 

RETURN 

C 

C  ENDIF 

C 

C 

c 

C22222222222222222222222222222222222222222222222222222222222222222222222 

C 

IF  (MAPTYPE  .EQ.  2)  THEN 
C  FIND  THETA  PROJECTIONS  ONTO  PLANE 
THETAN  *  DAT AN ( TRDI $T)/TRDI ST 
C 

XOUT  *  TANPHI  *  THETAN 
YOUT  ■  TTHETA  *  THETAN 
C 

RETURN 

ENDIF 

C 

C22222222222222222222222222222222222222222222222222222222222222222222222 

C 

C 

033333333333333333333333333333333333333333333333333333333333333333333333 

C 

C 


0034 


IF  (  MAPTYPE  .EQ.  3) THEN 
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C 

r 

FIND  TAN-THETA  PROJECTION  WITH  PRIMARY  DISTORTION 

0035 

l 

FIRST  =  TANPHI*TANPHI  +  TTHETA*TTHETA 

0036 

r 

DTERM  =  1  +  DFACTOR  *  FIRST 

0037 

t 

XOITT  =  TANPHI  *  DTERM 

0038 

r 

YOUT  -  TTHETA  *  DETERM 

0039 

V 

RETURN 

0G40 

ENDIF 

C 

0333333333333333333333333333333333333333333333333333333333333333333 

C 


V 

C 

0444444444444444444444444444444444444444444444444444444444444444444 

r 

0041 

IF  (MAPTYPE  .EQ.  4)  THEN 

c 

r 

FIND  SIN  THETA  PROJECTIONS  ONTO  PLANE 

0042 

t 

r 

SINN  =  DS IN ( DATAN ( TRDI ST ) ) /TRDIST 

0043 

XOUT  *  TANPHI  *  SINN 

0044 

r 

YOUT  =  TTHETA  *  SINN 

0045 

r 

RETURN 

0046 

t 

ENDIF 

0444444444444444444444444444444444444444444444444444444444444444444 

C 

c 

0047  RETURN 

0048  END 


INPUT  : 


OUTPUT  : 


uni  iiiwiiiiiiiiww.il  i.i|iiwiii)inii|||jjpgB—B 

{> 


NAVTRAEQUIPCEN  IH-332 

SUBROUTINE  FOV 


-  VIA  PLOTTING 
ROUTINES 


FIELD  OF  VIEW  #1  (FCVl) ,  FIELD  OF  VIEW  #2  (FCW2) ,  NUMBER  OF  POINTS 
ON  FO/'s  (NPT) ,  LINE  TYPE  FOR  FOVl  (Ll) ,  LINE  TYPE  FOR  F0V2  (L2) , 
DEGREES  TO  RADIANS  FACTOR  (CCNVR)  ,  MAPTYPE  (MAP) ,  PRIMARY  DISTORTION 
FACTOR  (DFACTOR) . 

TO  GRAPHIC  TEPMINAL  ONLY  -  CIRCIDS  INDICATING  FCVl  AND  FCV2. 
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TESTS  L'B.  FOR.  1 

C 

c 

c 

c 


C 

C 

★  * 

C 

*  FOV  * 

C 

*  * 

C 

C 

C 

C 

THIS  SUBROUTINE  GENERATES  POINTS  FOR  FOV  BOUNDARIES 

C 

CORRESPONDING  TO  F0V1  AND  FCV2,  FOP  TANGENT  MAPPING 

C 

MAP  =1,  FOR  IDEAL  R-THETA  MAPPING  MAP  =2,  FOR  P-THETA 

C 

MAPPING  WITH  PRIMARY  DISTORTION  MAP  =  3,  FOR  SIN  THETA 

c 

c 

MAPPING  MAP  =  4. 

c 

c 

r 

A  MAXIMUM  OF  75  POINTS  IS  ALLOWED 

L 

c 

THIS  ROUTINE  IS  CALLED  AFTER  THE  GRAPHIC  TERMINAL  IS 

c 

c 

INITIATED  BY  A  CALL  TO  CPLOT  (IE.  SCALING  IS  SET) 

0001 

c 

r 

SUBROUTINE  FOV  (FOV1,FOV2,NPT,L1,L2,CONVR,MAP,DFACTOR) 

0C02 

1/ 

r 

DIMENSION  XI (75) ,X2(75) ,Y1 (75) ,Y2(75) 

L 

c 

r 

FIND  RADII  OF  BOUNDARY  CIRCLES 

0003 

L 

IF  (  MAP  .EQ.  1  )THEN 

0004 

R1  »  TAN(  FCV1 *0 . 5  *  CONVR  ) 

0005 

R2  -  TAN(  F0V2*0.5  *  CONVR  ) 

0006 

r 

ENDIF 

0007 

L 

IF  (  MAP  .EQ.  L  )THEN 

c 

FIND  RADII  FOR  R-THETA  BOUNDARIES 

oooe 

R1  *  F0V1 *C . 5  *  CONVR 

0009 

R2  -  F0V2*0.5  *  CONVR 

0010 

r 

ENDIF 

0011 

l 

IF  (  MAP  .EQ.  3  )THEN 

c 

FIND  RADII  FOR  TANTHETA  W/PRIMARY  DISTORTION  BOUNDARIES 

0012 

R1  =  TAN(  F0V1 *0.5  *  CONVR  ) 

0013 

R2  *  TAN(  F0V2*0.5  *  CONVR  ) 

0014 

R1  ■  RT1  *  (  1  +  DFACTCR  *  RT1*RT1  ) 

0015 

R2  *  RT2  *  (  1  +  DFACTOR  *  RT2*RT2  ) 

0016 

ENDIF 

c 

IF  (MAP  .EQ.  4)  THEN 
R1  *  SIN(  FOV1  *  .5  *  CONVR) 


'  ( 


0017 

0018 
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0019 

0020 


0021 

0022 

0023 

0024 

0025 


0026 

0027 

0028 

0029 

0030 


0031 

0032 


0033 

0034 

0035 

0036 
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R2  =  SIN(  F0V2  *  .5  *  CONVR) 

ENOIF 

C 

C  FIND  RADIAN  INCREMENT  FOR  POINTS  ON  CIRCLE 
C 

RADINCR  =  2.  *  3.14159265349  /  (  NPT  -  1  ) 

C 

DO  100  N  =  l.NPT 
C 

ARGH  *  (N  -  1 )  *  RADINCR 
CSA  =  COS (ARGH) 

SNA  *  SIN(ARGH) 

C 

C  STORE  POINTS  ON  CIRCLE 

C 

XI (N)  =  R1  *  CSA 

Y1  N)  =  R1  *  SNA 

X2(N)  =  R2  *  CSA 

Y2(N)  =  R2  *  SNA 

C 

100  CONTINUE 
C 

C  DRAW  F0V1  WITH  LINETYPE  'Ll* 

C 

CALL  SETLINE(Ll) 

CALL  LINE  (XI, Y1 . NPT, 1, -1,0. ,1,0,0) 

C  DRAW  F0V2  WITH  LINETYPE  ' L2 ' 

C 

CALL  SETLIME(L2) 

r  CALL  LINE(X2, Y2, NPT, 1, -1,0., 1,0,0) 

RETURN 

END 
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<1-  ~ 
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0001 

0002 

0003 


0004 

0005 

0006 

0007 

00G8 

0009 

0010 

0011 
0012 
CGI  3 
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C 

C 

C 

C 

c  *  * 

C  *  TBLDIST  * 

C  *  * 

C 

C 

0  DISTORTION  PERCENTAGES  AR  FOUND  ACCORDING  TO  DEFINITION'S  FOUND 
C  IN  THE  IRE  STANDARDS  FOR  TELEVISION. 

C 

C 

C  THIS  SUBROUTINE  IS  INTENDED  TO  BE  USED  WITH  DISTORTION  ROUTINES 
C  THAT  UTILIZE  A  ONE  DIMENSIONAL  ARRAY  OF  POINTS  THAT  DEFINE 

C  A  SET  OF  POINTS  ON  A  PLANAR  SURFACE.  THE  ALGORITHM  DEPENDS  ON 

C  BEING  SUPPLIED  WITH  THE  'X'  AND  'Y'  POINTS  FOR  DISTORTION 
C  CALCULATIONS  AND  THE  'AC'  AND  'YC'  POINTS  FOR  COMPARISON. 

C 

C  NVERT  -  NUMBER  OF  POINTS  DOWN  THE  PLANE 

C  NHORIZ  -  NUMBER  OF  POINTS  ACROSS  THE  PLANE 

C 

C  CONSIDERING  A  RECTANGULAR  PLANE  ,  THE  POINTS  THAT  ARE  USED  ARE: 
C  THE  POINTS  AT  THE  FOUR  CORNERS,  AT  THE  MIDDLE  OF  THE  OUTER 

C  EDGES,  AND  AT  THE  CENTER  OF  THE  PLANE.  (  NINE  POINTS  ) 

C  THEREFORE  THE  VALUES  NVERT  AND  NHCRIZ  MUST  BE  ODD. 

C 

c 

c 

SUBROUTINE  TBLDIST  (X,Y, AC, YC, NVERT, NHORIZ) 

C 

INTEGER*4  1(9) 

DIMENSION  PCTCI ST ( 9 ) ,X(1),Y(1),XC(1),YC(1) 

C 

C  CHECK  NVERT  AND  NHORIZ  FOR  ODDNESS 
C 

PNHD2  =  NHORIZ  /  2.0 
NHD2  =  NHORIZ  /  2 
RNVD2  =  NVERT  /  2.0 
NVD2  =  NVERT  /  2 
C 

HN  =  RNHD2  -  NHD2 
VN  =  RNVD2  -  NVD2 
C 

ODD  «  0.5 
C 

IF(  (HN  .NE.  ODD)  .OR.  (VN  .HE.  ODD)  )  THEN 
TYPE  *  ,'  NUMBER  OF  POINTS  IS  NOT  ODD,  TBLDIST  EXITED' 
RETURN 
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0014 

C 

C 

r 

ENDIF 

CALCULATE  INDEX  VALUES  FOR  THE  NINE 

0015 

l 

Ki) 

a 

1 

0016 

1(2) 

= 

RNHD2  +0.5 

0017 

1(3) 

NHORIZ 

0018 

1(4) 

- 

NVD2  *  NHORIZ  +  1 

0019 

1(5) 

1(4)  +  1(2)  -  1 

0020 

1(6) 

= 

1(5)  +  1(2)  -  1 

0021 

1(7) 

= 

(NVERT  -  1)  *  NHORIZ  +  1 

0022 

1(8) 

= 

1(7)  +  1(2  -  1 

0023 

r 

1(9) 

s 

1(8)  +  1(2)  -  1 

V 

C 

r 

FIND  MAGNIFICATION  RATIO  ’ACONST' 

0024 

L 

PKT 

a 

Y ( I (2 ) )  -  Y(I(8) ) 

C025 

CPHT 

a 

YC(T(2) )  -  YC( I (8 ) ) 

0026 

ACONST 

=  PUT  /CPHT 

0027 

0028 

0029 

0030 

0031 

0032 


C 

C 

C 

C 


C 

C 

C 

C 

C 

C 


LOOP  FOR  CALCULATING  PERCENT  DISTORTION  FOR  THE  NINE  POINTS 


DO  100  N  *  1,9 

IN  =  I(N) 

TX  =  X ( IN) 

TY  =  Y(IN) 

TXC  *  XC(IN)  *  ACONST 
TYC  =  YC(IN)  *  ACONST 


FIND  RADIAL  DISTANCE  FROM  POINT  ON  VIEW  PLANE  TO  POINT 
ON  VIEW  WINDOW  ' GPE 1 . 


r 

GPE  =  SQRT((TXC-TX)*(TXC-TX)+(TYC- 

V 

c 

r 

CALCULATE  DISTORTION 

PCTDIST(N) 

1  *  GPE  /  PHT  *  100. 

19 

ni 

CONTINUE 

0036 

c 

WRITE (4,*) 

i,' DISTORTION  PERCENTAGES 

0037 

r 

WRITE(4,2) 

0038 

WRITE(4,3] 

* ( I ( L ) ,L=1 ,3) 

0039 

WRITE(4,1 

,{PCTDIST(J),J®1,3) 

0040 

WRITE  4,2 

0041 

WRITE  4,3] 

,(I(L),L*4,6) 

0042 

WRITE(4,1 ] 

, (PCTDIST(J) ,J-4,6) 

s 
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0043 

WRITE(4,2) 

0044 

WRITE(4,3  ,{l(L) ,LS7,9) 

0045 

c 

WRITE(4,1 )  , (PCTDIST ( J) ,Js7,S 

0046 

V 

p 

RETURN 

C047 

i 

F0RMAT(5X,F6.2,V,3X,F6.2,'% 

0048 

3 

FORMAT  4X.I3.6X, 13, 6X, 13) 

0049 

2 

P 

FORMAT(/) 

0050 

L 

END 
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